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Abstract 

We use rigorous non-equilibrium thermodynamic arguments to prove (i) the residual entropy of 
any system is bounded below by the experimentally (calorimetrically) determined absolute tem- 
perature entropy, which itself is bounded below by the entropy of the corresponding equilibrium 
(metastable supercooled liquid) state, and (ii) the instantaneous entropy cannot drop below that 
of the equilibrium state. The theorems follow from the second law and the existence of internal 
equilibrium and refer to the thermodynamic entropy. They go beyond the calorimetric observa- 
tions by Johari and Khouri [J. Chem. Phys. 134, 034515 (2011)] and others by extending them 
to all non-equilibrium systems regardless of how far they are from their equilibrium states. We 
also discuss the statistical interpretation of the thermodynamic entropy and show that the con- 
ventional Gibbs or Boltzmann interpretation gives the correct thermodynamic entropy even for a 
single sample regardless of the duration of measurements. 
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I. INTRODUCTION 



In a recent publication, Johari and Khourii used upper and lower bounds on the isobaric 
entropy S{Tq) of a glass as a function of the temperature Tq of the surrounding medium, 
obtained by using the measured heat capacity during cooling and heating to argue for 
the reality of the residual entropy Sn, as the latter has become a highly debated issue in 
the literature.-"— For a brief review of the history of the residual entropy and the current 
controversy, we refer the reader to Goldstein,-- Gutzow and Schmelzer,'i2 Nemilov,~ and 
the recent reviews^i^^ by us; see also below. The existence of a non-zero residual entropy 
{Sji > 0) is very common in Nature, and does not violate Nernst's postulate, as the latter 
is apphcable only to equilibrium states with a non- degenerate ground state; see Sect. 64 
in Landau and Lifshitz.^^^ Its existence was first theoretically demonstrated by Pauling and 
Tolman;— see also Tolman.— In addition, the existence of the residual entropy has been 
demonstrated rigorously for a very general spin model by Chow and Wu.— The residual 
entropy for glycerol was observed by Gibson and Giauque^^ and for ice by Giauque and 
Ashley. Pauling2^ provided the first numerical estimate for the residual entropy for ice, 
which was later improved by Nagle.— Nagle's numerical estimate has been recently verified 
by simulation.— The numerical simulation carried out by Bowles and Speedy^^ for glassy 
dimers also supports the existence of a residual entropy. Richet^- uses the Adam-Gibbs 
theory to justify the residual entropy. Thus, it appears that the support in favor of the 
residual entropy is quite strong. We wish to emphasize that what is customarily called the 
third law due to Nernst, according to which the entropy must vanish at absolute zero, is 
merely a postulate and not a strict theorem even in equilibrium.— i^i^ Indeed, many exactly 
solved statistical mechanical models show a non- zero entropy at absolute zero. However, as 
of yet, no experiment can be performed at absolute zero to demonstrate the residual entropy; 
in all cases, some sort of extrapolation is required. This point should not be forgotten in the 
following whenever we speak of measuring the residual entropy. In addition, we will speak 
of the equilibrium state associated with a non-equilibrium state. Depending on the context, 
the equilibrium state may represent a true equilibrium state such as a crystal or a (time- 
independent) metastable state such as the supercooled liquid. However, for the purpose 
of clarity, we will consider the supercooled liquid in the following, but the arguments are 
applicable to both cases. 
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A. Controversy and Its Current Status 



It is surprising to see this controversy persist in the current hterature even though it 
seemed resolved a long time ago;— we also note somewhat recent attempts.—*^ The con- 
troversial issue is the following: As the irreversibility does not allow for an exact evaluation 
of the entropy, is it possible for the entropy to decrease by an amount almost equal to S-r 
within the glass transition region so that the glass (see Glass2 in Fig. [1]) would have a van- 
ishing entropy at absolute zero? The entropy for Glass2 follows from the recent proposals by 
Kivelson and Reiss,- and by Gupta and Mauro;-'^'^^ see also Reiss.— Gupta and Mauro^ 
conclude that 

S < SscL (1) 

at and below the glass transition, in conformity with Glass2. They actually state the above 
inequality at the glass transition in Eq. (10) of the above paper but then take it to be also 
valid below the transition. Their conclusion is based on the fact that G > Gscl for the Gibbs 
free energies and the equality H = ifscL of the enthalpy at the glass transition; here S, G and 
H refer to the glass and ^scl, Gscl and ifscL refer to the supercooled hquid. The equality 
H = HscL is widely accepted as an experimental fact in the field.— The motivation for 
their proposal is their under standing^i^ii^i^^i^^ that the entropy of a single sample of glass 
at absolute zero must be zero as the microstate of the sample will not change no matter 
how long the glass is "observed." Thus, according to this view, the entropy is not merely 
less than that of Glassl, such as Glass3, it must be strictly less than or equal to ^scl- The 
alternative shown as GlassS 

•S* < ^Glassl > 'S'scL (2) 

is not consistent with this view.— 

According to Eq. ([1]), the entropy of a glass can only increase during relaxation as the 
glass strives to equihbrate at any temperature below the glass transition. What could be 
wrong with such a simple deduction, which does not require any sophisticated mathematics? 
Gutzow and Schmelzer— (see the discussion following their Eq. (39)) also come to the same 
conclusion. On the other hand, the common view of the glass is shown by Glassl in Fig. [U 
with the clear implication that the entropy of the glass is higher than that of the supercooled 
liquid. There is at present no consensus as is evident from the discussion reported in the 
proceedings,- and the general conclusion drawn by Goldstein^ that the abrupt entropy loss 
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FIG. 1: Schematic behavior of the entropy: equihbrated supercooled hquid (sohd curve) without 
any relaxation and three possible glasses (Glassl-dotted curve, Glass2-dashed curve, Glass3-dash- 
dotted curve) during vitrification as a function of the temperature Tq of the medium. Structures 
appear to freeze (over an extremely long period of time) at and below Tqq; see text. The transition 
region between Tog and Tqg over which the liquid turns into a glass has been exaggerated to high- 
light the point that the glass transition is not a sharp point. For Tq < Tog, the non-equilibrium state 
undergoes isothermal structural relaxation in time towards the supercooled liquid. For Glassl, the 
two vertical downward arrows show isothermal structural relaxation at two diff'erent temperatures, 
during which the entropy decreases. The same will also happen for Glass 3. For Glass2, the en- 
tropy increases during isothermal structural relaxation; see the upward arrow. The instantaneous 
temperature T{t) of the glass decreases towards the Tq during relaxation., so that the entropy is a 
function of T{t) during relaxation. The entropy of the supercooled liquid is shown to extrapolate 
to zero per our assumption, but that of Glassl to a non-zero value and of Glass2 to zero at absolute 
zero. The entropy of Glass 3 may or may not vanish at absolute. The possibility of an ideal glass 
transition, which does not affect our conclusion, will result in a singular form of the solid curve. 
The second law is used to support Glassl and rule out Glass2; see the text. 

comparable to S'r as proposed by Kivelson and Reiss^*^ will result in the construction of a 
perpetual motion machine of the second kind. Goldstein merely follows the consequence of 
the existence of a possible reversible connection between a glass and the supercooled liquid as 
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proposed by Kivelson and Reiss. In view of the simple observation in Eq. ([T]), one can argue 
that the general demonstration by Goldstein is either inapplicable to a glass or that the basic 
premise of a possible existence of a reversible connection between a glass and the supercooled 
liquid using Kivelson-Reiss construction must be invalid. Indeed, Gupta and Mauro^ii^ argue 
that the second law does not apply to the glass transition proposed by them as the loss 
of entropy in their formulation does not result in a latent heat, while Goldstein's general 
demonstration noted above depends crucially on this heat and, therefore, does not invalidate 
the proposal of a discontinuous loss by Gupta and Mauro although it does invalidate the 
proposal by Kivelson and Reiss. One can also argue that the special "electrode" construction 
proposed by Kivelson and Reiss, and used by Goldstein, requires microscopic information 
to confine the system into a unique basin, a configurational microstate, in that one needs to 
know precisely where each particle of the system is located within a certain small volume 
cell associated with its possible vibrations;—*^ see Gujrati^^ for elaboration on microstate 
measurements. Oppenheint^ has also raised somewhat of a similar objection. In addition, it 
is not obvious that thermodynamics, which is a description of a macrostate, can be applicable 
to a microstate such as the one obtained by the above special construction. 

Mauro et al^ also consider a glass formed by continuous cooling such as Glass2 in which 
the entropy is lost continuously. It should be noted that Goldstein^"^ does not explicitly 
discuss continuous cooling; however, in a private discussion, Goldstein has pointed out that 
his demonstration based on solubility also covers continuous cooling. The amount of loss 
in this case is not going to be close to S'r. Therefore, Goldstein's argument will not only 
invalidate Glass2 but also GlassS. 

Despite all these attempts, the situation remains confusing, in part because, to quote 
Goldstein^-: "As the residual entropies found by calorimetric measurements either equal 
or at least do not exceed the calculated value...., it is generally though not unanimously 
accepted that the residual entropies are real." Thus, in these cases, as already pointed 
out by Bestul and Chang,— one "...cannot demonstrate..." whether the residual entropy 
"...differs from zero." Of course, the attempts so far by those who believe in the residual 
entropy, ^''^li^'i^i^ have been only to demonstrate that the residual entropy is approximately 
equal to the calculated value, to be denoted here by 5'expt(^o); within a few percentages; 
however, as they do not deal directly with S'scl, they have not ruled out the inequality in 
Eq. ([T]). Even though no experiments have revealed a situation in which S'expt(^o) li^s below 
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S'scl(To) over the temperature range where the latter has been experimentally determined, 
the above calorimetric studies make no claim that 

5'cxpt(^o) > Sscl{To) (3) 

for all temperatures below the glass transition, even though many workers believe it to be 
intuitively true. In light of Eq. ([1]), however, t is possible that the above inequality turns out 
to be incorrect. The contribution of Goldstein^ only proves that an abrupt or a continuous 
loss of a part of the residual entropy {S < Sciassi) cannot be valid, but it neither rules out 
zero residual entropy nor the inequality in Eq. ([1]). It also says nothing about the above 
inequality. Is it possible that the demonstration of S* < ^scl itself is invalid? However, we 
cannot find any discussion of it in the current literature. Thus, it is not surprising that the 
controversy still persists. Goldstein^ has suggested that the entropy loss proposal " ...may be 
impossible to verify by any conceivable experiment,..." thus leaving the possibility that some 
theoretical approach, such as the one to be taken here, may resolve the issue. A theoretical 
approach will also allow us to investigate the inequalties in Eqs. ([UlS]) at absolute zero, 
where experiments cannot be performed. 

B. Clausius Limits 

Johari and Khourii have analyzed data from a large number of glasses to draw the con- 
clusion that a non-zero residual entropy is real in the cases they have analyzed. The idea 
behind their discussion is simple: the bounds, called the Clausius limits,- on the (thermody- 
namic) entropy reflect the contributions of irreversibility in the experimentally determined 
estimate SexptiTo) of the entropy S(To). They find that the two bounds are so tight in most 
cases (less than 2%) that one can neglect the corrections to the calorimetrically determined 
S'expt(To) due to irreversibility along cooling and heating paths, a conclusion arrived at by 
several others.-*^*^"— It appears that this approach was first used by Bestul and Chang^ 
and later by Sethna and coworkers.—. Johari and Khouri finally conclude that the exper- 
imental evidence of a non-zero residual entropy (see Glassl in Fig. [1]) is beyond reproach 
for the systems they have analyzed. The implication of their contribution and of several 
others mentioned above^'i^i'^^'^ is that there is no thermodynamic reason, such as the third 
law, for the residual entropy to be zero, a conclusion well known in theoretical physics.— 
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We should mention at this point that if there is ever any conflict between the second law^- 
and any other law in physics such as the zeroth or the third law for a macroscopic body, it 
is the second law that is believed to hold in all cases. This suggests that Glass2 does not 
materialize for the systems considered by Johari and Khouri.- 

C. Importance of Bounds 

We wish to prove in this work that Glass2 cannot materialize for any system. As far as 
Glass3 is concerned, we will establish that it may represent an approximation for Glassl, but 
cannot represent a real glass. Therefore, we will mostly consider Glassl and Glass2 in the 
following. To fulfill this goal, we have to go beyond the calorimetric evidencei'^'i^'i^ for the 
residual entropy. We treat the two glasses in Fig. [T]basically as two separate proposals for the 
general behavior of glasses. Therefore, they need to be demonstrated to be valid or invalid in 
all cases, without a single exception. To obtain such a general result, we must not make any 
assumptions or approximations; the latter would mean that the conclusion could not be valid 
in all cases; see below. We must also not rely upon heuristic or intuitive arguments as part of 
the proof. As the exact values of the entropy require detailed information about the system, 
it in not feasible to find the two entropies exactly for all systems. Thus, we will not be 
interested in system-specific knowledge; rather, we want mathematically sound conclusions 
about the entropy that will be valid for all glass formers without exception. Because the 
controversy is between two inequalities, any approximation would turn the strict inequality 
such as 5 < SscL into S^^^'^"'^ ^ 5'scl- In that case, we may have 5* > ^scLin some cases 
depending on how far S is from 5'^^^''°^. Unfortunately, various discussionsi'^i'i^ii^ in favor 
of the residual entropy and the specific critique by Goldstein^ using solubility either use 
approximations (such as replacing the inequality in Eq. below by an equality) or are 
phenomenological with a limited domain of validity. Indeed, Gutzow and Schmelzer— have 
come to the conclusion (see above Eq. (40) there) that in any real process, the entropy 
either remains constant or increases. The conclusion is in accordance with Eq. [1] Thus, the 
current status of the field leaves open the possibility that the residual entropy could vanish. 

It is abundantly clear from the above discussion that there is a need to look at the issue of 
residual entropy once again. In our opinion, the real issue is the inequality in Eq. ([T]), whose 
validity has not been scrutinized in the literature. If the bound is jusitified, it automatically 
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proves that S'r is not real. If, on the other hand, the bound is not justified, it does not 
automatically prove that S'r is real. We then have the additional task to prove its reality. 
As is customary, we treat the supercooled liquid as an equilibrium state, even though it not 
a true equilibrium state; see above. We proceed by following the strict second law inequality 
d[S > in Eq. (|H]) and use it along with the existence of internal equilibriunt^ (for Theorem 
|2]) to prove the following two theorems applicable to all non-equilibrium systems, regardless 
of how close or far they are from their equilibrium state. 

Theorem 1 The experimentally observed non-zero entropy at absolute zero in a vitrification 
process is a strict lower bound of the residual entropy of any system: 

SR = SiO)> S,,ptiO)> Sscl{0). (4) 

Theorem 2 Any drop of the glass entropy below that of the supercooled liquid (such as 
Glass2 in Fig. \^ is a violation of the second lawM Thus, 

S > SscL, (5) 

so that the entropy variation in time has a unique direction as shown by the downward 
arrows in Fig. [H 

Our conclusion in Eq. ([5]) is in contradiction with that in Eq. ([1]). Therefore, we need 
to understand the reason for this discrepancy. All experiments on or exact/approximate 
computations for non-equilibrium systems must obey the strict inequalities in Eqs. (jUlS]) 
without any exception. This is the meaning behind the usage of "... rigorous ..." in the title 
of the paper. The actual values of the entropy are not relevant for the aim of this work, which 
is to settle the controversy between Glassl and Glass2 under vitrification and the way their 
entropies relate to that of the equilibrated supercooled liquid. Because of the possibility that 
the systems may be far away from equilibrium, where the irreversible contributions may not 
be neglected, our results go beyond the previous calorimetric evidence.-'^'^'i^ The systems 
we are interested in include glasses and imperfect crystals as special cases. However, to be 
specific, we will only consider glasses below; the discussion is valid for all non-equilibrium 
systems. 

We hope that the proofs of the above two theorems settle the controversy between the 
two forms of glasses and about the residual entropy S'r. 
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We first consider the behavior of the thermodynamic entropy that appears in classical 
non-equilibrium thermodynamics. As is well known, this entropy is governed by the second 
law.— In the next section, we prove Theorem [TJ In the following section, we prove Theorem 
121 Both sections deal with the thermodynamic entropy. In Sect. IIV[ we explain how the 
thermodynamic residual entropy can be understood in terms probability of a microstate, 
and how even a single sample can give a highly reliable result for the residual entropy. The 
last section contains our conclusions. 

II. FORWARD ENTROPY BOUND DURING VITRIFICATION: THEOREM [1] 

The process we consider is carried out at some cooling rate as follows. The temperature 
of the medium is isobarically changed by some small but fixed ATq from the current value 
to the new value, and we wait for (not necessarily fixed) time Xobs at the new temperature to 
make an instantaneous measurement on the system before changing the temperature again. 
At some temperature Tog, see Fig. [H the relaxation time r^einx, which continuously increases 
as the temperature is lowered, becomes equal to Tobs- Just below Tog, the structures are 
not yet frozen; they "freeze" at a lower temperature Tqg (not too far from Tog) to form an 
amorphous solid with a viscosity close to 10^^ poise. This solid is identified as a glass. The 
location of both temperatures depends on the rate of cooling, i.e. on Tobs- Over the glass 
transition region between Tqq and Tog shown in Fig. [H the non- equilibrium liquid gradually 
turns from an equilibrium supercooled liquid at or above Tog into a glass at or below Tqg, a 
picture already known since Tammann.^^ Over this region, some dynamical properties such 
as the viscosity vary continuously but very rapidly. However, thermodynamic quantities 
such as the volume or the enthalpy change continuously but slowly. As the observation 
time Tobs is increased, the equilibrated supercooled liquid continues to lower temperatures 
before the appearance of Tog. In the hypothetical limit Tobs — )■ C)0, it is believed that the 
equilibrated supercooled liquid will continue to lower temperatures without any interruption, 
and is shown schematically by the solid blue curve in Fig. [H We overlook the possibility 
of the supercooled liquid ending in a spinodal.— It is commonly believed that this entropy 
will vanish at absolute zero (5'scl(0) = 0), as shown in the figure. As we are going to be 
interested in S'scl(^o) over (0, Tog), we must also acknowledge the possibility of an ideal glass 
transition in the system. If one believes in an ideal glass transition, then there would be a 
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singularity in S'scl(To) at some positive temperature Tk < Tqq, below which the system will 
turn into an ideal glass whose entropy will also vanish at absolute zero.— The possibility of 
an ideal glass transition, which has been discussed in a recent review elsewhere,— will not be 
discussed further in this work. All that will be relevant in our discussion here is the fact that 
the entropy vanishes in both situations (5'scl(0) = 0). However, it should be emphasized 
that the actual value of S'scl(O) has no relevance for the theorems. 

It is a common practice to think of the glass transition to occur at a point that lies 
between Tog and Tqq. Gupta and Mauro^ consider the glass transition to occur at Tqq 
to obtain the bound in Eq. ([T]). We will not make this assumption in this work except 
when we discuss their inequality later. We have drawn the two entropy curves (Glassl or 
Glass2) in Fig. [T] that emerge out of 5'scl(^o) for a given Tobs in such a way that Glassl 
has its entropy above (so that > 0) and Glass2 below (so that Sn = 0) that of the 
supercooled liquid. The entropy of Glassl (Glass2) approaches that of the equihbrated 
supercooled liquid entropy from above (below) during isothermal (fixed temperature of the 
medium) relaxation; see the two downward vertical arrows for Glassl. It is the approach 
to equilibrium that distinguishes the two glasses, Glassl and Glass2. Although Johari and 
Khouri do not mention in their conclusion, their analysis of tight bounds also shows that 
the entropy does not drop by an amount close to within the glass transition region for 
the systems studied by them. However, because of the involved approximation, it sheds no 
light on whether the glass entropy lies above or below the corresponding S'scl(7o). Thus, 
their work and many others leave open the possibility that Glass2 may materialize if the 
irreversibility is too large. This again shows why obtaining a bound is so important, even if 
we do not determine the actual entropy values. 

The concept of internal equilibriunt^'^ for a non-equilibrium system means that its 
instantaneous entropy is a state function of its instantaneous state variables like energy, 
volume etc. and any internal variables^'^"— used to specify its state. Their usage is also 
a common practice^"^ now-a-days for glasses. Employing internal equilibrium gives rise 
to an instantaneous Gibbs fundamental relation, see Eq. ( |T^ below, which determines its 
instantaneous temperature, pressure, etc. 

We now prove Theorem [1] Consider an isobaric process (we will not explicitly exhibit 
the pressure in this section) from some state A at temperature Tq in the supercooled liquid 
state which is still higher than Tog to the state Aq at absolute zero. The state Aq depends 
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on the path A— j-Aq along with Tq = 0, which is imphcit in the following. We have along 
A^Ao 

S{0) = S{To)+ Jd,S+ Jd,S, (6) 

A A 

where we have assumed that there is no latent heat in the vitrification process, and 
where^i^"— dS = deS + diS, each of which for a non-equilibrium system is path depen- 
dent. The component 

deSit) = -dQ{t)/To = CpdTo/To (7) 

represents the reversible entropy exchange with the medium in terms of the heat dQ{t) 
given out by the glass at time t to the medium whose temperature at that instant is Tq. The 
component d^S > represents the irreversible entropy generation in the irreversible process. 
In general, the irreversible term 

diS > (8) 

contains, in addition to the contribution from the irreversible heat transfer with the medium, 
contributions from all sorts of viscous dissipation going on within the system and normally 
require the use of internal variables.— i^""— Thus, Eq. (j6]) contains all possible sources of 
entropy variations.— i^"— This is easily proven by considering the system and the medium as 
an isolated system^^ and expressing the entropies as functions of the instantaneous values 
of the observables and internal variables. A discontinuous change in the entropy is ruled out 
from the continuity of the Gibbs free energy G and the enthalpy H in vitrification proved 
elsewhere.— Thus, we only consider a continuous change in the entropy as shown by the two 
glass curves in Fig. [H 

The equality in Eq. (jH]) holds for a reversible process, which we will no longer consider 
unless stated otherwise. The strict inequality d^S > occurs only for irreversible process 
such as in a glass. Since the second integral in Eq. (|6]) is always positive, and since the 
residual entropy S'r is, by definition, the entropy S{0) at absolute zero, we obtain the 
important result 

Sn = S{0) > 5expt(0) = S{To) + JCpdTo/To. (9) 

To 

This confirms the expectation noted above that the irreversibility during vitrification does 
not allow for the determination of the entropy exactly, because determining the second inte- 
gral in Eq. (|6]) is not trivial, especially if internal variables are not considered.— The for- 
ward inequality is due to the irreversible entropy generation from all possible sources^i^"— 
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that seems to not have been recognized by the proponents^'^'i^i'^i^^ of vanishing S'r. This 
strict forward inequahty clearly establishes that the residual entropy at absolute zero must be 
strictly larger than S'expt(O) in any non-equilibrium process. This proves the first inequality 
in Eq. (HD. 

We now prove the second inequality in Eq. We consider processes that occur when 
Tohs < Treiax(To). Let Q{t) = dQ(t)/dt bc the rate of net heat loss by the system. Then, for 
each temperature interval dTo < below Tog, we have 



'^obs 



\dQ\ =Cp\dTo\ = J 







Q 



TrclaxCTb) 

dt<\dQl^{To)^ J 





Q 



dt, Tn < Tn, 



where |(iQ|gq (To) > denotes the net heat loss by the system to come to equilibrium, i.e. 
become supercooled liquid during cooling at Tq. For Tq > Tog, we have dQ = dQeq(TQ). 
Thus, 

JCpdTo/To > J Cp,eq(iTo/To, 

To To 

where |(i(5|eq = Cp^eq \ dTo\. We thus conclude that 

So.pt{0) > ^scl(O); (10) 

the strict inequality is the result of the fact that glass is a non-equilibrium state. Otherwise, 
we will have 5'cxpt(0) > 5'scl(0) for any arbitary state. 
This proves Theorem [H 

The difference S'r— 5'expt(0) would be larger, more irreversible the process is. The quantity 
S'expt(O) can be determined calorimetrically by performing a cooling experiment. We take Tq 
to be the melting temperature Tqm, and uniquely determine the entropy of the supercooled 
liquid at Tqm by adding the entropy of melting to the crystal entropy S'cr(Tom) at Tqm- The 
latter is obtained in a unique manner by integration along a reversible path from Tq = to 
To = Tom: 

Tom 

'S'cr(Tom) = 'S'cr(O) + / Cp^cRdTo/To, 



here, 5'cr(0) is the entropy of the crystal at absolute zero, which is traditionally taken to 
be zero in accordance with the third law, and Cp^cr(To) is the isobaric heat capacity of 
the crystal. This then uniquely determines the entropy of the liquid to be used in the 
right hand side in Eq. ([9]). We will assume that S'cr(O) = 0. Thus, the experimental 
determination of 5'expt(0) is required to give the lower bound to the residual entropy in Eq. 
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(jl]). Experiment evidence for a non-zero value of 5'expt(0) is abundant as discussed by several 
authors various textbooks^iiS^ also discuss this issue. Goldstein^ gives a value of 
Sr ~ 15.1 J/K mol for o-terphenyl from the value of its entropy at Tq = 2 K. However, we 
have given a mathematical justification of S'expt(O) > in Eq. f lTU]) . The strict inequality 
proves immediately that the residual entropy cannot vanish for glasses, which justifies the 
curve Glassl in Fig. [H 

The inequality in Eq. (|9]) takes into account any amount of irreversibility during vitri- 
fication; it is no longer limited to only small contributions of the order of 2% considered 
by Johari and Khouri and by several others which makes our derivation very gen- 
eral. The relevance of the residual entropy has been discussed by several authors in the 
literature.-a^-i^^-i^^ 

By considering the state Aq above to be a state Aq of the glass in a medium at some 
arbitrary temperature Tq below Tog, we can get a generalization of Eq. (|9]): 

T' 

5(T^) > 5expt(T^) = Sin) + jcpdn/n. (n) 

To 

We again wish to remind the reader that all quantities depend on the path A— t-Aq, which 
we have not exhibited. By replacing Tq by the melting temperature Tqm and Tq by Tq, 
and adding the entropy ^(Tom) of the medium on both sides in the above inequality, and 
rearranging terms, we obtain (with 5'l(Tom) = 5'scl(Tom) for the liquid) 

sunu) + sinu) < Sin) + sinu) - / Cpdnin. (12) 

Tom 

where we have also included the equality for a reversible process. This provides us with an 
independent derivation of the inequality given by Setna and coworkers.— 

It is also clear from the derivation of Eq. ffTOj) that the inequality can be generalized to 
any temperature Tq < TQg with the result 

^-exptlTo) > 5'scl(To), (13) 

with S'expt(ro) 'S'scl(^o) as ^0 ng from below. Thus, S'expt(ro) appears in form similar 
to that of GlassS in Fig. [H except that the latter represents a possible glass entropy while 
the former represents the calorimetric approximation for Glassl. 

While we have only demonstrated the forward inequality, the excess S'p — »S'cxpt(0) can 
be computed in non-equilibrium thermodynamics,—!^"— which provides a clear prescription 
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for calculating the irreversible entropy generation. We do not do this here as we are only 
interested in general results, while the calculation of irreversible entropy generation will, of 
course, be system-dependent and will require detailed information. Gutzow and Scmelzer— 
provide such a procedure with a single internal variable but under the assumption of equal 
temperature and pressure for the glass and the medium. However, while they comment that 
diS > whose evaluation requires system-dependent properties, their main interest is to 
only show that it is negligible compared to deS. 

We have proved Theorem [1] by considering only the system without paying any attention 
to the medium but assuming the second law as is evident from Eq. (jH]). We have done 
this because the proponents of vanishing S'r normally consider the glassy state without ever 
bringing in the medium in the discussion. This does not mean that the conclusion would 
be any different had we brought the medium into our discussion. This is seen from the 
derivation of the ineqaulity in Eq. ( !T2|) from Eq. (ITT]) . We will find it convenient to consider 
the medium in the next section to overcome the objection^i^ that the glass does not obey 
the second law. 

III. ENTROPY AND ENTHALPY DURING RELAXATION 

We now turn to the inequality in Eq. ([T]) to see if it would ever be satisfied. To avoid 
directly discussing the relationship of the latent heat with the entropy loss, the possibility of 
entropy loss, or whether the second law applies to a glass, we change our mode of presentation 
and consider the system not by itself, but as a part of an isolated system in which the system 
is surrounded by an extremely large medium whose temperature Tq, pressure Po, etc. are not 
affected by what happens within the system. To prove Theorem [2l we consider the system to 
be not in equilibrium with the medium. All processes that go on within the medium occur 
at constant temperature, pressure, etc. Thus, there will not be any irreversible process 
going on within the medium. All irreversible processes will go on within the system. This 
simplification occurs because of the extremely large size of the medium and will be central 
in our discussion here. 

While some may doubt that the second law is not applicable to the glass at the glass 
transition, there cannot be any doubt that the second law applies to the isolated system. It 
is found^*^"^ that the instantaneous temperature T{t), pressure P{t), etc. of the system are 
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different from the corresponding quantities of the medium when the former is not in equihb- 
rium with the medium. All that is required is for the system to be in internal equilibrium^^ 
which is defined as a state in which the entropy has no explicit time-dependence; its time 
variation is due to time- dependent observables and internal variables. The existence of 
instantaneous T(t), P(t), etc. is a consequence of internal equilibrium, and is a general 
property of any system out of equilibrium even at high temperatures and is not restricted 
to glasses only. Therefore, T{t), P(t), etc. should not be confused with fictive temperature 
and pressure, etc. that are meaningful for glasses. This issue has been discussed earlier.— 
The Gibbs fundamental relation for the system when it is in internal equilibrium is given by 



where we have allowed for a single internal or structure variable ^ (t) for the sake of simplicity; 
for glasses, see Nemilov^ and Gutzow and Schmelzer— for the usage of internal variables. We 
have explicitly shown the time- dependence in the above equation to highlight the presence of 
relaxation in the system. The affinity A{t) is conjugate to the internal variable and vanishes 
when the system comes to equilibrium with the medium {Aq = 0). In that case, we also 
have T{t) Tq and P{t) Pq; see Eq. (fT5|). 

For a system out of equilibrium, the instantaneous entropy S{t) and volume V{t) seem to 
play the role^ of "internal variables," whose "affinities" are given by T(t) — Tq and P(t) — Pqj 
respectively. This fact is not common in the glass literature to the best of our knowledge. 
The temperature and pressure of the system are usually taken to be those of the medium, 
which is an approximation. For example, Schmelzer and Gutzow^ identify dcS = CpdT/T, 
see their Eq. (1), whereas it should be properly identified as in Eq. ([7]) with T replaced by 
Tq. They also identify the pressure in their Gibbs fundamental relation [see their Eq. (2)] 
as the external pressure. 

We now turn to prove Theorem |2i Let us rewrite the Gibbs fundamental relation as 



dE{t) = TodS{t) - PodV{t) + (T(t) - To)dS{t) - (P(t) - Po)dV{t) - A{t)d^{t), (15) 



in which each of the last three terms can be associated with an irreversible entropy 
generation.— For this, it is easier to take all but the first term on the right side to the 
other side of the equation. We thus note^^ that 



dE{t) = T{t)dS{t) - P{t)dV{t) - A{t)di{t) 



(14) 



(To - T{t))dS{t) > 0, (P(t) - Po)dV{t) > 0, A{t)d^{t) > 0, 



(16) 
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in accordance with the second law. The equahties above and below occur only for reversible 
processes. As we are only interested in irreversible processes in non-equilibrium systems, 
the inequalities above and below become strict inequalities, which cannot be violated in any 
real process. Thus, as before, we will exploit these strict inequalities to derive a bound on 
the rate of entropy variation. 

We extend the derivation given earlier^^ to include the internal variable to obtain as the 
statement of the second law:— 

dt \T{t) To) dt \T{t) To) dt T{t) dt - ' ^ ' 
each term in the first equation must be non-negative. In a vitrification process, in which 
the energy decreases with time, we must, therefore, have 

T{t) > To 

during any relaxation (at a fixed temperature and pressure of the medium) so that T{t) 
approaches Tq from above [T{t) — ?■ T^] as the relaxation ceases and the equilibrium is 
achieved. It now follows from Eq. (fT6l) that during vitrification 

dS{t)/dt < 0; (18) 

the equality occurring only when equilibrium with the medium has been achieved. The 
above inequality gives the bound on the rate that we are interested in. At the end of the 
relaxation 

S{To,Po,t) — )■ S^Q.i^{To, Po); 

the plus sign is to indicate that the glass entropy reaches 5'scl(7o, Pq) from above. 

We have shown To,Po in S(To,Po,t) = S{T(t), P(t), A(t)) to emphasize that the result 
is general during any relaxation. In the derivation, which only considers the behavior of 
So{t) of the isolated system, no assumption about the nature of irreversibility such as any 
loss of ergodicity in the system, inapplicability of the second law to the system, possibility 
of any chaotic behavior, chemical reaction, etc. is made. The only assumption that has 
been made is that it is possible to define the instantaneous temperature and pressure for the 
system. We have also made no assumption that S{t) lies above (Glassl, Glass3) or below 
(Glass2) the entropy Sscl(To) of the equilibrated supercooled liquid; see Fig. [H Being a 
general result, it should be valid for any real glass. This now gives a way to decide which 
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of the glasses in Fig. [T]is in accordance with the above conclusion. Above Tog, the system 
is always in equilibrium with the medium so its temperature is the same as Tq. Below Tog, 
when the system is not in equilibrium with the medium, then T{t) > Tq as long as there 
is no equilibrium. The entropy of Glassl and Glass3 approach Ssch{To) of the equilibrated 
supercooled liquid from above during any isothermal relaxation, which is consistent with 
Eq. ( !T8|) . As the entropy is a unique function of the path, the two glasses must correspond 
to different histories. Therefore, from now on, we will not consider GlassS anymore. 

For Glass2 in Fig. [H the entropy actually drops below that of the supercooled liquid 
by some amount;"i2^i^ the amount does not even have to be comparable to 5*^. In this 
situation, the entropy must approach Sscl(Tq) of the supercooled liquid from below during 
relaxation; see the upward arrow in Fig. [T] This will result in the increase of the entropy 
during relaxation, which violates Eq. f|T8|) . Thus, Glass2 cannot be rationalized. For the 
same reason, the conclusion of Gutzow and Schmelzer— of increase in entropy cannot be 
rationalized. 

We are now ready to investigate what could be technically wrong with the entropy loss 
proposal. This will also settle whether the glass must obey the second law or not at the glass 
transition. From the behavior of dSQ{t)/dt in Eq. f|T7j) . we can immediately identify^'^'^^ the 
rates for the entropy of the system and the medium, respectively, 

dSjt) _ 1 dE{t) Pit) dVjt) Ajt) djit) 

dt ~ T{t) dt ^ T{t) dt ^ T{t) dt 
dSjt) _ 1 dEjt) _ Pq dVjt) 

dt Tq dt To dt 

While the entropy change of the medium has no irreversible contribution as noted earlier, 
the irreversible entropy change diS{t) of the system is given by the three terms in Eq. ([T7j) . 
each of which must be non-negative. Writing dS{t) = deS(t) + diS(t), we find that 

deS{t) _ dS{t) _ 1 dE{t) PodV{t) 
dt dt Tq dt Tq dt 

d,S{t) dSQjt) ^ ^ 

dt dt ~ 

In general, dcS{t)/dt can have either sign. In a cooling process, we have dE{t)/dt < 0. 
Moreover, we normally have dV{t)/dt < 0. Thus, deS{t)/dt < 0. However, its sign is not 
relevant for our discussion. On the other hand, d[S{t)/dt > 0, which follows from the second 
law apphed to the isolated system, for which there is no dispute. With these consequences 
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of the second law, we can now evaluate the merit of the entropy loss proposal.-*^'^'^'^ Let 
us consider vitrification. The change in the entropy dcS{t) < is due to exchanges with 
the medium. This is part of the entropy change that will occur even in a reversible process. 
The contribution from the entropy loss,-'^'^^!^'^ which we denote by dSioss{t), is due to 
the vitrification process. Vitrification in the entropy loss view represents changes occurring 
within the system. As the contribution from every internal process must be included in 
diS{t), and as each such contribution must be non-negative, there is no way to justify the 
negative contribution dSioss{t) due to vitrification^'^'i^i^'^ without violating the second law 
for the isolated system. This proves Theorem [2J 

This is the conclusion obtained by Goldstein;^ we have just provided a direct proof of 
his conclusion. There is also no merit to the suggestion of Gupta and Mauro^'i^i'^^ that the 
glass does not obey the second law. 

It is easy to see that the discussion above can be easily extended to include other ob- 
servables and internal variables without affecting Theorem O The above discussion was 
also not restricted to a constant pressure of the medium. Indeed, the discussion above has 
been very general. The only restriction was the extremely large size of the medium, which 
is easily satisfied in experiments. For the general case, the Gibbs free energy is given by 
G{t) = Hit) - ToS{t). At fixed Tq, we have^ 



dGit) _ dH{t) dS{t) 



from which it follows that 



dt dt 

dH{t) 



dt 



dt 
dS{t) 



(19) 



dt 



We now turn to considering isobaric vitrification during which the system is not in equilib- 
rium with the medium. We will assume that the system is always very close to mechanical 
equilibrium so that its pressure is equal to Pq; however, there is normally no thermal equi- 
librium so that the instantaneous temperature T{t) of the system is different from Tq. We 
will now show that the above general conclusion remains unaltered. From Eq. f lT^ . we have 



dH{t) 



dt 



>T{t) 



dS{t) 



dt 



(20) 



where we have used A(t)d^(t)/dt > 0. The last bound is tighter than the previous bound 
and reduces to the equality obtained earlier— in the absence of any internal variable ^. 
In any case, the enthalpy of the glass is not constant in time at Tqq even for an isobaric 
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vitrification if we accept tliat tlie entropy undergoes relaxation. We tlius conclude that there 
is no justification in assuming AH{Tqg) = 0. 

We now turn to the inequality in Eq. ([T]) to inquire what may be wrong with its trivial 
justification offered by its authors. We consider Tq < Tog- The enthalpy H{Tq) of the glass 
relaxes towards its equilibrium value Hscl(Tq) of the supercooled liquid from above during 
isobaric vitrification, contrary to the assumption by Gupta and Mauro.^ Acknowledging this 
immediately leads to 

^og(5'(Tog) — 5'scl(^og)) < H{Tog) — HscLiToc). 

As the right side is strictly positive, there is no justification in concluding the inequality in 
Eq. (II]). It is most certainly possible to satisfy the above inequality and also have 

'S'(Tog) > 'S'sclI^og)- 

This is consistent with strict inequality in Eq. (fT8|) . The equality AH(Tq) = is only valid 
at Tog, where the difference AS(Tq) = 0, so that S(Tog) = SschiTog). 

The isothermal relaxation that occurs during glassy vitrification originates from the ten- 
dency of the glass to come to thermal equilibrium during which its temperature T{t) ap- 
proaches To in time. The relaxation process results in the lowering of the corresponding 
Gibbs free energy^ in time, as expected; this results in not only lowering the enthalpy dur- 
ing vitrification, as observed experimentally, but also of the entropy S{t) during relaxation, 
as shown for Glassl in Fig. [H 

IV. STATISTICAL INTERPRETATION OF THERMODYNAMIC 5r 

Now that we have established the reality of the residual entropy by considering the 
thermodynamic entropy in classical thermodynamics, we wish to discuss its possible statis- 
tical interpretation. We recall that the thermodynamic entropy cannot be given a unique 
value; all that we can discuss in thermodynamics is the change in it. It does not even have 
to be non-negative, as is evident from the entropy of an ideal gas at very low temperature. 
However, as we have defined S(To, t) with respect to S'cr(O) in this work, and we have taken 
'S'cr(O) = 0, the thermodynamic entropy has a unique value. Therefore, any attempt to 
provide a statistical interpretation must result in an agreement with the numerical value 
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of the above thermodynamic entropy, a point also made by Sethna and coworkers.— The 
issue of the statistical interpretation has been discussed elsewhere by us.— As a glass is 
a frozen structure (over an extremely long period of time), we index each frozen structure, 
which represents a microstate, at absolute zero hj i = 1,2, ■■■ ,Wg- Each microstate is 
characterized by the set of observables and internal variables. All glasses formed under 
identical macroscopic conditions will be in one of these microstates at Tq = 0; let pi > 
denote the probability that a glass will be in the microstate i. The instantaneous statistical 
entropy is an average quantityi^ii^ (we set fee = l):S(t) = — {Imp) = '^^pi{—\npi). If the 
glass formation occurs under an unbiased condition, all microstates will be equally probable 
{ Pi = p = I/Wq for all i) so that the residual entropy is = InlV^Q- The necessary 
(but not sufficient) condition for this is that the system be macroscopically large. The suf- 
ficient condition requires the system to be in internal equilibrium, though not necessarily 
in equilibrium with the medium, as discussed elsewhere.— In this case, the contribution 
— \np = InWc to the entropy from any of the Wq microstates is the same so that 

S = - \np = \nWG; (21) 

there is no need to sum over all microstates since YliPi ~ ^- J^^^ microstate will give 
us the correct entropy. Similarly, just one microstate will give a highly reliable value of any 
thermodynamic quantity. As the equiprobable condition will be overwhelmingly satisfied 
for a system in internal equilibrium, just one glass sample will give us a highly reliable 
thermodynamics. The discussion can be applied to a system at any temperature, provided 
it is in internal equilibrium. There is no need to carry out an ensemble or time average. This 
is what makes classical thermodynamics such a robust and highly reproducible endeavor, 
a result quite well known in equilibrium. We have just extended it to systems in internal 
equilibrium. 

Just because a glass sample at absolute zero is in a single microstate does not mean 
that its statistical entropy is zero. Such a statistical entropy interpretation does not agree 
with the thermodynamic entropy which, as discussed above, is known to yield a non-zero 
residual entropy in many cases. Therefore, such an interpretation, the one taken by Kivel- 
son and Reiss^ and by Gupta and Mauro,- must be considered physically irrelevant. The 
correct statistical formulation of the entropy is given by —\np, as shown in Eq. ( 12T|) . 
and merely reflects the sample probability of a glass prepared under identical macroscopic 
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conditions.— Any sample will give the same statistical entropy equal to the thermodynamic 
residual entropy so that the statistical and thermodynamic interpretations of the entropy 
are equivalent. 

The discussion should not imply that the ensemble average is no longer meaningful just 
because we are dealing with a single sample. The way to see it most clearly is to imagine 
dividing the system into a macroscopically large number of quasi-independent subsystems, 
each of which itself is macroscopically large. The subsystems are identical in size. Let l 
denote a microstate of one of these subsystems, and > its probability. Then, it follows 
from their quasi- independence that p = HPt! product is over all subsystems. It is now 

k 

easy to see that 

where the sum is over all subsystems; each subsystem can be considered as representing a 
member of the ensemble. All these members are considered at the same instant. 

Recently, Goldstein-^ has also discussed the relevance of a single microstate for the average 
at any temperature, but he does not discuss the time required for the system to come to 
internal equilibrium or to equilibrium after adding the enthalpy increment to it. Thus, it is 
by no means clear that a single microstate during this time will really represent the average 
properties of the system. His requirement of "...overwhelming majority of microstates..." 
must refer to the equilibrium state of minimum Gibbs free energy at given Tq,Pq, but 
this is not relevant for a glass. As the enthalpy increment is added, it takes a while for 
this enthalpy to distribute itself throughout the system from the boundary. During this 
interval, the system becomes inhomogeneous and will not even be in internal equilibrium, 
let alone in equilibrium. In this interval, we can still treat each subsystem discussed above in 
internal equilibrium. But then we are dealing with an ensemble average over inhomogeneous 
subsystems. We believe that his conclusion would be valid only after the system has come 
into internal equilibrium, though not necessarily in equilibrium, as demonstrated above. For 
this, we must replace his above requirement by equiprobability of microstates. 

V. CONCLUSIONS 

The current work was motivated by the confusion about the residual entropy and about 
the behavior of the entropy during relaxation that exists in the literature, as discussed in 
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Sect. II A[ There are opposing views not only about the residual entropy but also about the 
impact of Goldstein's observation for the latter. The use of calorimetric data so far has been 
to demonstrate that the irreversibility during a glass transition is minimal so that Scxpt{To) is 
not different from the actual thermodynamic entropy S{Tq) of the glass. But the calorimetric 
evidence does not reveal how S'(ro) or S'oxpt(^o) relate to 5'scl(^o)- Therefore, how S(Tq) 
approaches Sscl(To) remains unsettled; there are competing views in the literature. To 
clarify the situation, we have considered the role of irreversible entropy generation during 
isobaric vitrification and prove Theorems [1] and [2] that are valid regardless of how far the 
system is out of equilibrium, as long as it is in internal equilibrium; the latter is required 
to define the instantaneous temperature, pressure, affinity, etc. The theorems are very 
general and are not restricted by the " amount" of irreversibility. Theorem [T] shows that the 
calorimetrically measured absolute zero entropy forms a strict lower bound to the residual 
entropy. The former is shown to be positive under the assumption S'scl(O) = 0; otherwise, 
we have the strict inequality in Eq. ( ITOj) for a glass. In general, we have S'cxpt(O) > 5'scl(0). 
It then follows that the residual entropy has to be larger than 5'expt(0). Theorem [2] shows 
that the instantaneous entropy S{To,t) must always be higher than or at most equal to the 
entropy 5'scl(^o) of the equilibrated supercooled entropy, which invalidates the inequality in 
Eq. ([1]). During isothermal relaxation, the entropy must decrease towards 5'scl(^o) in time. 
All physical systems must follow the two inequalities in Eqs. (jUlS]) without any exception. 
We have demonstrated that the entropy loss proposal violates the second law and have put 
the original observation of Goldstein on firmer grounds. We have not only justified but 
also strengthened the calorimetric evidence of the residual entropy by establishing > 
5'expt(0) > Sscl(O). 

The theorems follow from considering the thermodynamic entropy that appears in the 
second law for an isolated system. Thus, any attempt to provide a statistical version of 
entropy must satisfy these two consequences. We have shown that the conventional statistical 
entropy formulation is consistent with the thermodynamic notion of the residual entropy 
and that the equiprobability requirement explains why a single sample is sufficient to give a 
highly reliable thermodynamics of the system even when the latter is not be in equilibrium 
with the medium. 
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